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Abstract. In previous attempts to perform seismic modelling of pulsating 
subdwarf-B stars, various mode identification techniques are used with uncer- 
tain results. 

We investigated a method so far neglected in sdB stars, but very successful 
for Main Sequence pulsators, that is, mode identification from the line-profile 
variations caused by stellar pulsation. 

We report the calculation of time-resolved synthetic spectra for sdB stars 
pulsating with various combinations of pulsation modes; these calculations were 
carried out over appropriate ranges of effective temperature, surface gravity 
and helium abundances. Preliminary tests using these synthetic line-profile 
variations demonstrated their potential for mode identification by comparison 
with observation. 
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1. INTRODUCTION 

Subdwarf B (sdB) stars are believed to be low-mass (~ O.SMq) core helium burn- 
ing objects belonging to the Extreme Horizontal Branch (EHB) (Heber 1986). 
With thin and mostly inert hydrogen-rich residual envelopes, sdB stars remain hot 
(20000 < Teff < 40000) and compact (5 < logg < 7) throughout their EHB hfetime 
(Saffer et al. 1994); they eventually evolve towards the white dwarf cooling se- 
quence without experiencing the Asymptotic Giant Branch and Planetary Nebula 
phases of stellar evolution. While binary population synthesis calculations by Han 
et al. (2002, 2003) successfully demonstrate the formation of sdB stars through 
several possible channels, resulting models require further comparison with obser- 
vation. As sdB stars are believed to be responsible for the ultraviolet-upturn seen 
in the energy distributions of elliptical galaxies and spiral galaxy bulges (Yi et al. 
1997), understanding sdB star formation may provide an important diagnostic for 
studying galaxy evolution and formation. 

The discovery (Kilkenny et al. 1997, Green et al. 2003) that some sdB stars are 
nonradial pulsators means asteroseismology could be used to discern their internal 
structure and so constrain evolution models. But before any seismic modelling 
can be reliably attempted one must identify pulsation modes which are excited. 
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Assuming radial and 
non-radial pulsations can 
be modelled by spheri- 
cal harmonics Ynim, mode 
identification tries to as- 
sign values to the spherical 
wavenumbers n, I and m: 
n is related to the number 
of nodes of the radial dis- 
placement, I is the number 
of nodal lines on the stellar 
surface and m is the num- 
ber of such lines passing 
through the rotation axis of 
the star, as shown in Fig. 1. 
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Fig. 1. 1 = 3 pulsation mode with m from 3 (left) 
to (right) and i from 30° (top) to 90° (bottom) 
(GONG /NSO/ AURA /NSF) . 



2. THE PROBLEM OF MODE IDENTIFICATION IN SDB STARS 

Mode identification in subdwarf-B stars is far from being trivial. Several methods 
which prove efficient for mode identification in pulsating Main Sequence stars have 
limitations in the case of sdB stars. Nevertheless, these mode identification tech- 
niques, however uncertain, are used in meritorious attempts at seismic modelling 
of the most promising pulsating sdB stars. 

2.1. Period fitting 

The method entails fitting an observed period spectrum of a pulsating star to a 
theoretical period spectrum and is straightforward in principle. However, when 
considering pulsating sdB stars there is a major difficulty to overcome (Brassard 
et al. 2001): so far, most sdB stars (except for PG 1605-1-072 with more than 
50 periods (Koen et al. 1998), and KPD 1930+2752 with at least 44 periods 
(Billeres et al. 2000)) have a rather sparse period spectrum, possibly because 
lower amplitude frequencies are not yet observed. This makes it very difficult to 
match objectively observed periods with periods computed from a model, as the 
latter are much more numerous. 

Despite the intrinsic difficulty of period fitting in sdB stars. Brassard et al. 
(2001) use this technique to model the pulsating sdB star PG 0014-1-067 and ten- 
tatively identify 23 distinct pulsation modes. They search a four-parameter space 
(Tcff, log g, the mass of the star M, and the mass of its H-rich envelope Monv) for 
a model whose theoretical period spectrum could account for the observed periods 
in PG 0014-1-067, using a "merit" function that is minimum for the best fit(s). 
In addition to the four parameters, they also deduce the rotation velocity of the 
star and other properties (see their Table 6). They claim a good agreement with 
previous spectroscopic determination of Tcff and log g, but we note that in order 
to discriminate between equally good fits (several minima of their merit function) 
they use spectroscopic constraints. More recently, the same method is applied by 
Charpinet et al. (2005) to PG 1219-1-534. 

2.2. Amplitude ratios 

Yet another method in use for mode identification in pulsating Main Sequence 
stars relies on multicolor photometry. This method, known as the amplitude 
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ratio method (Hcyndcrickx et al. 1994), relies on the fact that for non-radially 
pulsating stars the photometric pulsation amplitude depends on wavelength and 
on the spherical wavenumber I of the pulsation mode. Ramachandran et al. (2004) 
investigate this method for non-radially pulsating sdB stars and draw the following 
conclusions: 

• for EC 14026 stars, only the modes of high I (3 and 4) should be easy to disen- 
tangle, provided their amplitudes are large enough to observe. On the other 
hand, low degree modes are very difficult to distinguish, and this problem 
becomes worse as the period shortens; 

• for "Betsy" stars, modes of spherical degree I = 1 may be difficult to distin- 
guish from modes with Z > 1 in the case of longer periods. 

Using ULTRACAM multicolor photometry, Jeffery et al. (2004, 2005) use the am- 
plitude ratio method to identify pulsation modes of three sdB stars, HS 0039+4302, 
KPD 2109+4401 and PG 0014+067. Their results are provisional, as only statisti- 
cal errors are taken into account. 

3. LINE-PROFILE VARIATIONS IN PULSATING SDB STARS 

We saw in the previous section that current mode identification methods in sdB 
stars need improvement. However none of these techniques use spectroscopy to 
perform mode identification, although it has proved useful when modelling pulsat- 
ing Main Sequence stars (e.g. Acrts et al. 1992, Briquet & Acrts 2003). Our claim 
is that spectroscopy, and the pioneering detailed computation of line-profile varia- 
tions (Ipv) we perform, could provide a more reliable mode identification method. 
Indeed Ipv enable a complete reconstruction of the non-radial oscillations of a star, 
provided that a time-series of spectra are observed with a high enough S/N ratio, 
temporal, and wavelength resolution. 

3.1. Computation of synthetic line-profile variations 

The method used for our computation of theoretical line-profile variations is quite 
straightforward. Wc combined throe computer programs, BRUCE (Townsend 
1997), SYNSPEC (Hubeny & Lanz 2000) and KYLIE (Townsend 1997), so as 
to compute a series of time-resolved spectra exhibiting line-profile variations. 

BRUCE first divides the stellar surface into a large number of very small sur- 
face elements that arc determined by the usual polarand azimuthal coordinate 
angles 9 and (p and a step-size A9 and A(p. BRUCE then perturbs the stellar 
surface, assuming the pulsation is linear and adiabatic, with selected combina- 
tions of pulsation modes at selected time steps, giving for each surface element its 
temperature, surface gravity, orientation, projected radial velocity and projected 
area. 

It is important to realize that, because of the pulsation, different points on the 
stellar surface not only have different radial velocity, but also different tempera- 
ture, log g and orientation, meaning that the contributions of the different surface 
elements to the line-profile have different amplitudes and energy distributions. The 
line-profile variations do not only come from the Doppler shift in wavelength, but 
also from the complex temperature and log g behavior on the stellar surface. 

The emergent monochromatic flux in wavelength interval AA, at some wave- 
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Fig. 2. Trial computation of line-profile variations for the short-period pul- 
sator PG 1605-1-072. Only the six highest-amplitude pulsation modes (out of more 
than 50) were taken into account in this computation in order to reduce computing 
time. The line-profiles have been shifted along the ordinate axis (normalized flux) 

for clarity. 

length A and time t is then given by 



where I\{t)ij is the specific intensity emergent from the tile at 6i and ipj in wave- 
length interval AA, at wavelength A and time i, in the direction of the observer, 
as computed by our modified version of SYNSPEC. The quantity 



accounts for the time-dependent orientation of each tile, n{t)ij being the normal to 
the tile at some time t, and s being the direction of propagation of the radiation. 

This is then repeated for the desired number of time-steps using KYLIE, there- 
fore producing series of time-resolved observer-directed spectra. 

Examples of line-profile variations computed with our codes are given for dif- 
ferent modes in Figs. 2 and 3. 

3.2. Observation of line-profile variations in pulsating sdB stars 

Observations of line-profile variations in sdB stars have been made by Telting & 

Ostcnsen (2004) who observed PG 1325+101 with the 2.5-m NOT telescope. Their 
observations showed strong evidence for a dominant radial pulsation mode, but 
most importantly drew attention to the potential of spectroscopic mode identifi- 
cation. Similar observations were made of the bright sdB star Balloon 090100001 
(Telting & Ostensen 2005). 
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Fig. 3. A time series of lA resolution synthetic spectrum ratios for PB 8783 

(short-period pulsator) and EC 21324-1346 (long-period pulsator) in left and right 
frames respectively. Second and subsequent synthetic spectra in each time series 
are divided by the first, the elapsed time being given in minutes. Time-series 
calculations were based on assumed pulsation modes (Z = 3, m = 0; / = 2, m = 
0; / = 1, m = 0; Z = 2, TO = 1; Z = 3, m = 0) and velocity amplitudes from Jeffery & 
Pollacco (2000) for PB 8783 and on modes (i = 1, m = 1; i = 2, m = 1; Z = 3, m = 
0) and amplitudes to match photometric variations of EC 21324— 1346. 

4. THE MOMENT METHOD AND ITS APPLICATION TO SDB STARS 
J^.l. The moment method 

Mode identification from line-profile variations could in principle be carried out 

by a direct comparison of observed and synthetic spectra but we propose to adapt 
to pulsating sdB stars the so-called moment method. It was first introduced by 
Balona (1986a,b, 1987), and later refined and applied to Main Sequence pulsators 
by Aerts et al. (1992), Aerts (1996) and Briquet & Aerts (2003). 

In this method, an observed line-profile is replaced by its first few moments, and 
their time dependence studied to identify the pulsation modes. The nth moment 
of a line-profile p{v) (where v is the line of sight velocity component corresponding 
to the displacement from the laboratory wavelength) is be defined as 



< >= 



(2) 



All the information contained in the line-profile can be reconstructed from the 
entire series of moments of nth order, but Briquet & Aerts (2003) showed that 
most of the information is contained in the first three moments: 



• the first moment < u > is the centroid of the line-profile in a reference frame 
with origin at the stellar center; 
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• the second moment < > is rclat 

• the third moment < v'^ > relates 
Telting 2003). 

Once the moments of the ob- 
served hne-profilc have been computed, 
their variations are compared with the 
time-dependence of theoretical moments 
(Aerts 1992). These depend on the 
wavenumbers I and m of the pulsation 
mode, but also on its amplitude Vp, on 
the inclination i of the star's axis of 
rotation to the observer, on its rota- 
tion velocity va and on the width cr of 
the intrinsic broadening profile, allowing 
a complete reconstruction of non-radial 
pulsation. One then selects the most 
likely set of parameters {I, m, Vp, i, vq, a) 
by minimizing a so-called discriminant 
(Briquet & Aerts 2003), which is a func- 
tion that takes into account the quality 
of the fit between the theoretical and ob- 
served moments. 

In Fig. 4, wc show the three normal- 
ized moments for a synthetic radial mode 
we computed. Even without computing 
the discriminant, the phase-dependence 
of the observed moments can provide 
some useful information. For instance 



to the variance of the line-profile; 
the skewness of the profile (see also 




Phase 

Fig. 4. First, second and third mo- 
ments for the synthetic line profile of a 
radial mode computed using our codes. 



• the pcak-to-peak amplitude of the first moment gives an idea about the 

overall velocity range due to the oscillation; 

• if < t;^ > can be described by a single sine function with twice the frequency 
of the first moment it means that m = 0, while if the mode is sectoral (|m| = 
I) < v"^ > behaves sinusoidally with the same frequency as the frequency of 
the first moment. 



4-2. Application to pulsating sdB stars 

The major difficulty to overcome in order to apply the moment method to pul- 
sating sdB stars comes from its current formulation (Briquet & Aerts 2003) in 
which the moment method neglects the variation of the specific intensity (i.e. 
51\{0i, ipj) = 0) due to changes in T^g and log 5 during the pulsation. This can be 
a good approximation when dealing with spectral lines that are not sensitive to 
temperature variations (such as silicon lines in jS Ccphei stars (Dupret et al. 2002) 
and (De Ridder et al. 2002)) but in the case of sdB stars, temperature effects 
should not be neglected. 

We tested the moment method in its current formulation on Ipv of for var- 
ious pulsation modes, but the identification was not always satisfactory. However 
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it was extremely difficult to choose a suitable wavelength range over which the mo- 
ments were to be computed. We are investigating a systematic way to make that 
choice, as well as whether the use of the narrower and much weaker helium lines 
might improve the mode identification. Another possibility is to take temperature 
effects into account using of higher-order moments, as Balona (1987) suggests. 

5. FUTURE PROSPECTS 

Our development of spectroscopic methods for pulsation mode identification in 
subdwarf-B stars was stimulated by the acquisition of observational data. A start 
on this was made in August 2005 using the South African Astronomical Observa- 
tory 1.9-m telescope to obtain time-resolved spectroscopy of "Betsy" stars (Schoe- 
naers & Lynas-Gray 2006). 

Furthermore, Kawaler & Hostler (2005) show internal rotation to be expected 
in sdB stars, the rate depending on distance from the stellar centre. Mode splitting 
due internal rotation cannot be identified using the usual formula applicable in 
the case of surface rotation and may have led to incorrect mode identification in 
previous works. Special attention therefore needs to be given to the modelling of 
this effect in pulsating sdB stars, even though the moment method shouldn't be 
affected much by differential rotation, contrarily to period fitting. 

However, mode identification in pulsating sdB stars is a real challenge with the 
methods discussed above. Other methods might prove worth investigating. 

The so-called IPS method aims at comparing the observed amplitude and phase 
variations across the line-profile with the theoretical variations (Schrijvers et al. 
1997 and Telting & Schrijvers 1997a, b). It was recently used by De Cat et al. 
(2005) to identify pulsation modes in monoperiodic SPB (Slowly Pulsating B) 
stars, together with the method of photometric amplitude ratios and the moment 
method (with moments up to the sixth order). 

Yet another mode identification method, Doppler Imaging, should be inves- 
tigated. The method is based on the idea that line-profile variations of rapidly 
rotating stars are a sort of "doppler image" of the stellar surface within the line- 
profile. This method has been challenged because of its apparent lack of physical 
foundation, but Kochukhov (2004) claims a reliable recovery of the surface pul- 
sation velocity structures can be achieved for all types of pulsation geometries 
accompanied by significant line-profile variations. It remains to be seen whether 
the line-profile variations in sdB stars are large enough. 

To summarize, this paper presents several options for meeting the challenge 
presented by sdB star pulsation mode identification by spectroscopic methods. 
All appear to be worthy of further investigation and progress made to date is 
briefly summarized. 
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